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Abstract 

We study the inverse scattering problem for electric potentials and magnetic fields 
in M. d , d > 3, that are asymptotic sums of homogeneous terms at infinity. The main 
result is that all these terms can be uniquely reconstructed from the singularities 
in the forward direction of the scattering amplitude at some positive energy. 

1 Introduction 



1.1 Background. There are different formulations of the inverse scattering problem for 
the multidimensional Schrodinger operator. The most complete and difficult problem is to 
find a one-to-one correspondence between potentials on M. d and the scattering data. For 
the solution of this problem see the paper by L. D. Faddeev [TU] and also and [33]. A 
review and further references can be found in [22] and j2H]. Here the difficult part is to find 
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necessary and sufficient conditions on the scattering amplitude f(u, uj\ A), u, to G E> d 1 , A > 0, 
such that the corresponding potential is local. 

A simpler problem of uniqueness and reconstruction of the potential from the scattering 
data has different settings. One of them is uniqueness and reconstruction from the high- 
energy limit of the scattering amplitude which was considered in jS] and [3] (see also the 
books [3122] and the papers [3 12])- In [5] it is used that the first Born approximation is 
essentially the Fourier transform of the electric potential. In the time- dependent method 
of 12] the X-ray or the Radon transforms of the electric and the magnetic potentials are 
reconstructed from the high-energy limit. 

A similar problem of uniqueness and reconstruction of the electric and magnetic po- 
tentials from the scattering amplitude at a fixed positive energy was considered in [21] for 
potentials of compact support, and in [2H El EH IH2] for potentials decaying exponentially 
at infinity. On the contrary, for general short-range potentials the scattering matrix at a 
fixed positive energy does not determine uniquely the potential. Indeed, in JB] examples -in 
three dimensions- are given of non-trivial radial oscillating potentials with decay as |x|~ 3 / 2 at 
infinity such that the corresponding scattering amplitude is identically zero at some positive 
energy. Moreover, in dimension two there are examples ^1] of potentials with a regular 
decay as \x\~ 2 at infinity that have zero scattering amplitude at some positive energy. Nev- 
ertheless, if two potentials have the same scattering amplitude at some positive energy and 
they coincide outside of some ball they necessary coincide everywhere |34j . 

Actually, the same problem appears in different settings. Thus, it is proven in ^3 E] 
that the scattering matrix at a fixed positive energy uniquely determines an exponentially 
decreasing perturbation of a stratified media. As another example, we mention that the scat- 
tering matrix at a fixed quasi-energy uniquely determines (see [3^]) time-periodic potentials 
that decay exponentially at spatial infinity. 

1.2 Setting the problem. We consider the related (at least in spirit) problem of the 
unique reconstruction of the asymptotics of the electric and magnetic potentials if only the 
singularity of the scattering amplitude f(u,uj; A) on the diagonal v = uo is known at some 
fixed energy A > 0. This problem was already studied in ^TJEElEi]; where the methods 
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of microlocal analysis were heavily used. We mention also a related paper [23] where using 
the method of [7j the asymptotics of the electric potential is reconstructed if the scattering 
amplitude is known on some (perhaps arbitrarily small) interval of positive energies. 

We show that this problem admits a quite elementary solution. It consists of two steps 
and both of them are in some sense well known. Let us explain them for the reconstruction 
of the leading terms of the electric V and magnetic A potentials. The first is the small 
angle approximation to the scattering amplitude f(v,uj; A) (see, e.g., [20J, §127) given by the 
formula 

/0, u; A) w e -^-i)/4 fc (d-D/2 (27r) -( rf -i)/2 f e ik(y, u - u) R ^ w; A; v) ^ k = A y 2j 

(1.1) 

where V = (V, A), 11^ is the hyperplane orthogonal to u, 

/oo 
(y(y + tu)-2k(w,A(y + tu)))dt (1.2) 
-oo 

and v — > u). The second is the classical inversion formula for the Radon transform. Actually, 
for any d we use always the two dimensional Radon transform. 

1.3 Scattering amplitude. Let us recall now the definition of the scattering amplitude. 
If the electric, V(x), and the magnetic, A(x), potentials are real- valued functions satisfying 
the assumption 

\V(x)\ + \A(x)\ + |divA(x)| < C(l + \x\)- p (1.3) 
where p > (d+ l)/2, then the corresponding Schrodinger equation 

(zV + A(x)) 2 ip + V(x)ifj = X = k 2 >0, (1.4) 

has a solution with the asymptotics 

tfj{x, u, A) = exp (ik(u, x)) + f\x\~ {d ~ 1)/2 exp (ik\x\) + o^x^ 1 ^ 2 ) 

as \x\ — > oo. Moreover, such a solution is unique. The coefficient f(v,uj; A) depends on the 
incident direction u of the incoming plane wave exp(ik(uj, x)), its energy A and the direction 
v = of observation of the outgoing spherical wave \x\^ t - d ^ 1 ^ 2 exp(zA;|x|). The function 
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f(u,oj; A) is known as the scattering amplitude. In terms of the scattering amplitude, the 
scattering matrix S(X) is defined by the formula 

(S(X)u)(u) = u(oj) + ie <r(*-3)/4 A (d-i)/4 (2 7r )-( rf - 1 )/ 2 f f( U) oj; X)u{u)dw. (1.5) 

A discussion of these results can be found, for example, in jS7]. On the contrary, the 
scattering amplitude can be formally defined by equation (|1.5|) via the scattering matrix 
(see subs. 2.1, for its time-dependent definition) for all short-range potentials when p > 1 in 
(JO}. Moreover, if the assumption 

\d a V{x)\ + \d a A{x)\ < C a (l + \x\)- p -W, p>l, (1.6) 

is valid for all a, then, as shown in jJJ, f(u,oj; X) is a C°°-function of u, oj G S^ -1 away from 
the diagonal v = oj. 

1.4 Main results. The diagonal singularity of the scattering amplitude contains all 
information about the behaviour of the potentials V(x) and A(x) as \x\ — > oo. This allows 
one to recover their asymptotics at infinity. 

Indeed, let us consider first the leading terms of the asymptotics of V(x) and A(x) 
when formulae (jl.lj) and (jl.2J) can be used. Given the leading singularity of the scattering 
amplitude on the diagonal, we can recover from (jl.lj) the function R(y,uj] X) as the inverse 
Fourier transform of this singularity in the hyperplane 11^. If A = 0, then, using (jl.2j) . we 
can reconstruct the asymptotics of V(x) by the inversion formula for the Radon transform. 
In the general case, we take into account that the integrals 

/oo 
V(y + tu)dt, uje§ d -\ yeU^, (1.7) 
-oo 

and 

/oo 
(u,A(y + tu))dt, ojeS d ~ l , yGlL, (1.8) 
-oo 

are, respectively, even and odd functions of oj G S^ 1 , which determine both these integrals 
separately. This allows us to reconstruct, as before, the electric potential V, but it is of 
course impossible to find the magnetic potential A from the integral (jl.8|) . Nevertheless, we 
show that this integral determines the magnetic field F = curl A. One cannot of course do 



4 



better because the scattering matrix, as well as the integral (jl.8|) . are invariant with respect 
to gauge transformations (see subs. 2.2). 

Assuming that the potentials are asymptotically homogeneous, we actually use, for all 
d > 3, the inversion of the Radon transform in two-dimensional planes. This is important 
for the reconstruction of the magnetic field. If d = 2, then the integral (jl.7|) can be equal 
to zero for all y ^ for a non-trivial homogeneous function V (see Counter-example 4.4). 
Therefore, in dimension 2 one cannot hope to recover the asymptotics of a potential from 
the forward singularity of the scattering amplitude. 

Our approach extends to the complete asympotic expansion. Assume that an electric 
potential 

oo 

V{x)^Y, V &) (°) 

3=1 

and a magnetic field 

oo 

F(x)~^2F j (x) (1.10) 

3=1 

are asymptotic sums of homogeneous functions (see subs. 2.3, for precise definitions) of 
orders — pf and — rj m ' ) , respectively, where 1 < p^ < pf* < ■ ■ ■ and 2 < < < ■ ■ ■. 
Suppose that all singularities in the forward direction of the scattering amplitude at some 
positive energy are known. We recover all functions Vj(x) and Fj(x) in a recursive way 
using at every step the inversion of the two-dimensional Radon transform. In particular, our 
results imply that all asymptotic coefficients Vj(x) and Fj(x) are uniquely determined by 
these scattering data. 

Compared to papers [T7J [TH] we do not assume that the degrees of homogeneity pj = 
1 +j, = 2+j. On the contrary, we recover these numbers from the singularity of the 
scattering amplitude. In particular, the assumption p^ > 2, > 3 imposed in ITS] 
turns out to be unnecessary. Another difference of our work is that we show that both sets 
of functions Vj and Fj can be obtained from the singularities of the scattering amplitude at 
only one energy (not at two as required in [TT1 ITH]). 

We proceed from a complete description of the diagonal singularity of the scattering am- 
plitude. Such explicit formulae were obtained in the papers IHO] and jHHl IHH] • These formulae 
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generalize (jl.2j) but formula (jl.lj) remains the same. This shows that actually it is natural 
to regard the scattering matrix as a pseudo-differential rather than an integral operator. 
This point of view was probably first explicitly adopted in [3] where the asymptotics of the 
eigenvalues of the scattering matrix was found. 

Under some minor additional a priori assumptions our results combined with those of [31] 
show that the electric potential (not only its asymptotic expansion at infinity) is uniquely 
determined by the scattering amplitude at some energy (see Theorem 4.6). This result im- 
plies that the uniqueness theorem does not require the super-power decay of the potential at 
infinity. Of course, it does not contradict the examples of [S] because we exclude oscillations 
of the potential at infinity and those of JT] because we assume that d > 2. 

The description of the diagonal singularities of the scattering amplitude is given in Sec- 
tion 3. The (inverse) problem of the reconstruction of the asymptotic expansions of V and 
F from these singularities is considered in Section 4. The main result is formulated as 
Theorem 4.2. 

Finally, we note that our approach extends naturally to long-range potentials. 

The inverse scattering for the Schrodinger equation has many important applications in 
atomic, molecular and nuclear physics. Also, the inverse scattering problem at a fixed energy 
for the Schrodinger equation is actually the same as the inverse scattering problem at a fixed 
frequency for the wave equation with a variable speed. This problem has many important 
applications in applied science and in engineering. For the applications of inverse scattering 
for the Schrodinger equation and the wave equation see (27j . 

2 Basic Notions 

2.1 Scattering theory. Here we recall some basic definitions of scattering theory. We 
consider the Schrodinger operator 

H = (iV + A(x)f + V(x) (2.1) 

with electric V(x) and magnetic A(x) = (A^(x), . . . , A^ d \x)) potentials in the space L 2 (M d ) 
where d > 2. Under the condition that V, A as well as the divergence divA of A are real and 
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bounded functions the Hamiltonian H is well defined and self-adjoint on the Sobolev class 
H 2 (IR d ). Let us denote by Hq = —A the "free" Hamiltonian corresponding to the case V = 
and A = 0. Under assumption (jl.3|) where p > 1, the positive spectrum of the operator H 
is absolutely continuous (and its negative spectrum consists of eigenvalues) and the wave 
operators 

W± = s- lim e itH e~ itHo (2.2) 

exist and enjoy the intertwining property HW± = W±H (see, for example, [19| or 28 ). It 
follows that the scattering operator 

S = W- (2.3) 

commutes with Hq and is unitary. 

We are now able to define the scattering matrix. Let be the unit sphere in M. d , 
R+ = (0, oo) and let L 2 (R+, L 2 (S d_1 )) be the L 2 -space of functions defined on R+ with 
values in L 2 (E> d ~ l ). Define the unitary mapping 

T : L 2 (R d ) — * L 2 (R + ,L 2 (S d - 1 )) 

by the equation 

{Tu){oo-\) = 2- 1 / 2 A (d - 2)/4 n(A 1/2 ^), where u{£) = (2n)~ d / 2 [ e - lM u(x)dx 

is the Fourier transform of u and the spectral parameter A plays the role of the energy of a 
quantum particle. Then (J-'Hqu)(X) = X(J 7 u)(X) and 

(JFS«)(A) = S , (A)(J r w)(A). 

The unitary operator S(X) : L 2 (E> d ^ 1 ) — > L 2 (S a!_1 ) is known as the scattering matrix at 
energy A. If assumption is valid for all a, then, as shown in pQ, S(X) is an integral 
operator with C°°-kernel s(u, u; A) away from the diagonal v = oj. Moreover, if p > (d+ l)/2 
in (jl.3|) . then S(X) is related to the scattering amplitude introduced in Section 1 by formula 
()1.5|) . For these results see, for example, [37]. 

2.2 Magnetic fields. It follows from (|2.2|) that the operators S and 5(A) are invariant 
under short-range gauge transformations, A \— > A — V0, where G C°°(lR d ) and 9 a 0(x) = 
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0(|a;| _e_ ' a ') for all a and some g > as \x\ — >• oo. Therefore it is natural to associate the 
scattering matrix 5(A) directly to the magnetic field F(x). However, since the definition of 
5(A) was given in terms of the magnetic potential A(x), we have to discuss here the relation 
between magnetic fields and potentials. Below we use freely three-dimensional notation for 
the curl and the divergence keeping in mind that in the general case A is a 1-form and F is 
a 2-form. By definition, F(x) = cmLA(x), that is in terms of components 

F^\x) = diA^\x) - djA {i \x). (2.4) 

On the contrary, the magnetic potential can be reconstructed from the magnetic field F(x) 
such that div F(x) = only up to arbitrary gauge transformations. It is not convenient 
to work in the standard transversal gauge (x,A tr (x)) = since even for magnetic fields of 
compact support the potential A tr (x) decays only as (a;]" 1 at infinity. 

Therefore we use the procedure suggested in jlU] of construction of the short-range mag- 
netic potential for an arbitrary magnetic field satisfying the condition 

\d a F(x)\ <C a (l + \x\)- r - lal , r>2, Va. (2.5) 

Let us introduce the auxiliary potentials 



noo d poo d 

/ sVF^^ijrfs, A®(x)=- sVF (,j) (ra)i,rfs. (2.6) 

Jl 3=1 J ° 3=1 



Note that is a homogeneous function of order —1, cmlA^x) = for x ^ and 
A tr = A reg + Aoq. Next we define the function U(x) for x ^ as a curvilinear integral 

U(x)= j {A^yldy) (2.7) 

taken between some fixed point Xq ^ and a variable point x. It is required that (jL ^ xo ,x, 
so that, in view of the Stokes theorem, for d > 3 the function U(x) does not depend on the 
choice of a contour T XOjX and grad[/(x) = A^x). Let us now choose an arbitrary function 
i] G C°°(lR d ) such that rj(x) = in a neighbourhood of zero, rj(x) = 1 for, say, \x\ > 1 and 
set 

A(x) = A reg (x) + (1 - riix^A^x) - U(x)gr&dr)(x). (2.8) 
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Then curlA(z) = F(x), A G C°°(M. d ) if F G C°°(R d ) and A(x) = A reg (x) for |x| > 2. 
Therefore it follows from definition ()2.6|) that under assumption f)2.5j) satisfies estimates 
(HUH with p = r - 1. 

Below we always associate to a magnetic field F the magnetic potential A by formulae 
()2.6|) - (|2.8|) and then construct the scattering matrix S(X) in terms of the Schrodinger 
operator ()2.1|) . If another short-range potential A satisfies the relation curlv4(x) = F(x), 
then necessarily the scattering matrices corresponding to potentials A and A coincide. This 
allows us to speak about the scattering matrix 5(A) corresponding to the magnetic field F. 

2.3 Asymptotic expansions. Let us introduce some notation. We denote by S~ p = 
S- p (R d ) the set of C°°(R d \ {0})-functions f(x) such that d a f(x) = 0(\x\~ p ~^) as |x| -> oo 
for all a. Then, the standard Hormander class is defined as, <S- p (lR d ) := C°°(R d ) n<S~ p (R d ). 
An important example of functions from the class S~ p are homogeneous functions / G 
C°°(R d \ {0}) of order -p such that f(tx) = t~ p f(x) for all x G M. d , i^O, and t > 0. 

Let functions fj G S^ pj where pj — > oo (but the condition p^- < p J+1 is not required). 
The notation 



oo 



(2.9) 

means that, for any N, the remainder 



AT 

f — ^^fjE<S~ p where p= min pj. (2-10) 

j= i J - N+1 

In particular, if the sum (J2.9)) consists of a finite number iV of terms, then the inclusion 
(|2.1(J|) should be satisfied for all p. A function / G C°° is determined by its asymptotic 
expansion (J2.9|) up to a term from the Schwarz class S = S~°°. 

It follows from formulae (J2.4j) and (J2.6j) - (J2.8|) that the magnetic field F(x) is homo- 
geneous of order — r' m ' < —2 if and only if the magnetic potential A(x) is homogeneous of 
order — p( m ) = —f( m ) + 1 < —1. Therefore A{x) is an asymptotic sum 

oo 

A(x)~J2M*) (2- 11 ) 
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of homogeneous functions Aj(x) of orders —Pj if and only if the representation (jl.lUj) holds 
for F(x) with homogeneous functions Fj(x) of orders — . Here the coefficients Fjix) and 
Aj(x) are related by formulae (|2.4|) and ()2.6|) - ()2.8|) : in particular, = rj — 1. 

Adding terms which are equal to zero, we can assume, without loss of generality, that the 
numbers p^- and A- in ()1.9|) and (|1.1U|) are related by the equality p^- = — 1. Anyway, 
by the reconstruction procedure equality of some terms to zero can never be excluded. Set 
V(x) = (V(x),A(x)). Then expansions (jl.9|) and (jl.lUJl are equivalent to the expansion 

oo 

V(*)^X>i(*) ( 2 - 12 ) 

i=i 

in homogeneous functions Vj(x) = (Vj(x), Aj(x)) of order — pj where pj = p^ = pf\ 

3 The Direct Problem 

3.1 Singularities of the scattering amplitude. Following j3HJE2il, we give in this section 
a complete description of the diagonal singularity of the integral kernel s(u,u;X) of the 
scattering matrix S(X) for a fixed energy A > 0. The dimension d > 2 is arbitrary. Our goal 
is to construct, for all N, an explicit function such that the difference 

s - s w e C pW (L 2 (S d - 1 ) x L 2 (S d - 1 )) where lim p(N) = oo. (3.1) 

N^oo 

Since s is a C°°-function away from the diagonal v — u, it suffices to construct in a 
neighbourhood of the diagonal only. This construction is given in terms of special approxi- 
mate solutions of the Schrodinger equation ([1.4)1 . Let us set 

u { f\x,0 = e^+^^b^^O, e = leil e (3.2) 

where the phase 

/■CO 

$±(x,£) = T &A(x±ti))dt (3.3) 
Jo 

and 

v 

n=0 
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The coefficients b± are denned here by the recurrent formulae b± (x, = 1 and 

POO 

&i n+1) (*, i) —T f ( ±\x ± t£, i) dt, (3.5) 
Jo 

where 

= 2i(A- V<S>±,Vb { ± } ) - Ab^ 

+ (|V$±| 2 -2(A,V$±) + V +\A\ 2 + idwA-iA<S>±)b { ± ) . (3.6) 

One can arrive to these expressions for the coefficients 6± ^ by plugging ()3.2)1 and (j3. 3)1 into 
the Schrodinger equation ()1.4j) where A = |£| 2 . This yields the transport equation for the 
function h*±\x,£). Seeking it in the form ()3.4)1 and equating coefficients at the same powers 
of (2i\£\)~ n } we obtain recurrent equations for the functions b±\ They are solved by formulae 
()3.5|) and ljH.fi)) . This is a standard construction of approximate solutions of the Schrodinger 
equation going back in the context of the present paper to ?5J. It is important that the 
remainder arising when function (j3.2j) is plugged into the Schrodinger equation satisfies, for 
larger N, better and better estimates as |x| — > oo everywhere except the forward direction 
x = £ for the function u+(x,£) and except the back direction x = — £ for the function 
u_ (x, £). At the same time the estimates of the remainder improve as the energy A — > oo. 

A complete description of the diagonal singularity of the function s(z/, uj; A) is given in 
the following theorem obtained in j3~%j I3T?]. 

THEOREM 3.1. Let estimates (ji.6| hold for all a. Let ojq G be an arbitrary point, 
be the plane orthogonal to ujq and Q = Q(u>o, 5) C § n_1 be determined by the condition 
(u, u } > 5 > 0. Set x = uj z + y, ye U^, 

h^(x,0=e^^b^(x,a (3.7) 



& {N) (y, u,u; A) = 2 1 (u + uj,u )h + (y,ku)h-(y,kuj) 



+ (2ik) 1 h + (y,kv){d z h-){y,kxj)- (d z h + )(y,kv)h-(y,kxj) 

- 2i{A(y),u )h+{y, kv) h- (y, fco;)] , h± = h { ± ] , k = A 1/2 , (3.8) 

and define 

s {n \u,uj;X) = (2ix)- d+l k d - 1 [ e ik{v >"- v) ^ N \y,v,u-\)dy. (3.9) 



Then relation (jff.i)) holds. 



11 



It can be also shown that the C^^-norm of the kernel (|3.1|) is 0(\~ q ( N ^) as A — > oo 
where q(N) — *> oo for N — > oo. Thus, all singularities of s(i/, cu; A) both for high energies and 
in smoothness are described by the explicit formulae (|3.8j) . ()3.9|) . 

Note that the amplitude (|3.8|) is invariant under a gauge transformation A i— > A — V0, 
/i± i— > e - ^ /i±, where G C 00 is an arbitrary function. 

Formula ()3.9|) gives the singular part s^ N \u, lu; A) of the kernel s(u,u;X) for u,u E. 
£1(uq,5). Since cj G § d_1 is arbitrary, this determines the function s( N '(v, u>; A) for all 
v,uj& E> d -\ 

3.2 Symbol of the scattering matrix. According to (J3.9j) it is natural to regard 
the scattering matrix S(X) as a pseudodifferential operator (on the unit sphere) determined 
by its amplitude. For our purposes it is convenient to reformulate Theorem 3.1 in terms of 
asymptotic series. First, replace e 1 **^'^ in (J3.7|) by its Taylor expansion. Then, rearrange the 
representation (J3.8|) collecting together terms of the same power with respect to V = (V, A). 
Thus, the amplitude a(y, is, uj; A) of the pseudodifferential operator 5(A) admits the expansion 
into the asymptotic series 

oo 

a(y, u, a «(^> ^ w ' A )' 

n=0 

where a (z/, u; A) = 2~ 1 (v + u,u>o), ai depends linearly on V, a 2 depends quadratically on 
V, etc. It follows from formulae ([3.5)1 . ()3.6|) that a„ G l S _ ^ _1 - )n for all n and, moreover, 



2ika l (y, u,u; A) = 2~ 1 (z/ + u, uj ) j \V(y + tu) + V(y - tu) 

- 2k{v,A{y + tv)) -2k{u,A(y-tuj))jdt, (3.10) 

up to a term from the class S~ p . 

Applications to the inverse problem require that the pseudodifferential operator 5(A) be 
determined by its symbol. Let us recall briefly the standard procedure (see, e.g., [211 EI]) 
of the passage in local coordinates from the amplitude to the corresponding (right) symbol. 
Let us consider the orthogonal projection p of Q on the hyperplane ILj and set 

a(y,u;\)~ J^ar 1 ^)-^ a(y,t( V ),uj; X) , y G LT^, (3.11) 

a 

12 



where the inverse mapping t = p~ l : p(£l) — > is given by the formula t(rj) = (77, (1 — M 2 ) 1 ^ 2 )- 
Note that we have taken here into account the factor — k in the phase in (|3.9|) which differs 
our definition from the standard terminology. Then 

s{v, to; A) = {2ir)~ d+1 k dr - 1 I e~ ik ^a{y, to; A) dy (3.12) 

so that the scattering matrix S(X) can be regarded as a pseudodifferential operator on S d_1 
with right symbol a(y,u; A). This leads to the following conclusion. 

THEOREM 3.2. Let estimates (ji.6|) hold for all a. Then the scattering matrix S(\) 
is a pseudodifferential operator on the unit sphere. Its right symbol a(y,u;X) admits the 
expansion into the asymptotic sum 

00 

a(y,uj;X) ~l + J2a n (y,uj;X), wG§" y G IT,, (3.13) 

n=l 

where a\ depends linearly on V, 02 depends quadratically on V, etc. Explicit expressions for 

all functions a n are obtained by putting formulae (jj.&j) and (jff.iij) together. In particular, 
we have that 

a n G S- (p - 1)n (3.14) 

and, moreover, 

&i := a 1 - R G S~ p , (3.15) 
where the function R is defined by formula . Thus, 

a - 1 - R G S~ p+1 - £ , e = min{p - 1, 1}, (3.16) 

so that R(y, uj\ A) is the principal symbol of the pseudodifferential operator S(\) — I . 

It is important that, for all n, the values of a n (y,u; A), y G IT W , y 7^ 0, are determined 
only by the values of V(x) and its derivatives on the line x = to + ty, t G M, which does not 
pass through the origin. 

Note that the symbol a(y, ui; A) can be recovered from the kernel s(v, 00; A) by the inversion 
of the Fourier transform (|3.12|) : 

a(y,u;X)= [ e ik ^s(t(r,),u; A) 7 fo) dr,, y G 11^, (3.17) 

13 



where 7 G C™(M. d x ) is an arbitrary function such that 7(17) = 1 in some neighbourhood of 
zero and 7(77) = for, say, \rj\ > 1/2. 

Let a be the function constructed in Theorem 3.2. Let us denote by T the mapping that 
sends V into the function a — 1. Of course it is defined up to a symbol from the class S~°°. 
Thus, we put 

T(y,u;\',V) = a(y,u;\)-l. (3.18) 
Moreover, we distinguish the linear part R of T and set 

Q(y, uj; A; V) = T(y, u; A; V) - R(y, 00; A; V), (3.19) 

so that 



Q(y, u; A; V) = a ± (y, u; A) + ^ a n(Vi ^ - ) 0- 



n>2 



We need the following simple property of the mapping Q. 



en 



PROPOSITION 3.3. Suppose that V® G <S^ ( , j = 1,2, w/iere p( 2 ) > pW > 1. Th 
Q(V« + V< 2 >) - Q(V«) G l S- p(2)+1 - £ , e = min{p« - 1, 1} > 0. 

Indeed, it follows from (J3.15j) that 

5l (V« + V( 2 )) - a!(V«) = Si(V( 2 )) G 5^ (2) . 

To consider the functionals a n (V) for n > 2, let us temporarily write them as a n (V) = 
a n (V, V, . . . , V) where a n (V^, V^, . . . , V^) is linear with respect to its arguments V*- 1 ), 
V( 2 ), . . . , VW. If V&'> G for all j = 1, ... , n, then 

a n (vW,..,v( n ))er (1) --" ( " )+n 

This property generalizes ([3.14)1 . Therefore, 

a2 ( V « + V ( 2 )) - a 2 (V«) = a 2 (V« V< 2 >) + a 2 (V( 2 \ V«) + a 2 (V^) 

and the right-hand side belongs to the class S~ p 1 ~ p2 +2 . The terms a n (V) for n > 2 can be 
considered quite similarly. 



14 



3.3 The case of asymptotically homogeneous potentials. Let us specify Theo- 
rem 3.2 for the perturbations considered in this paper. The following result is obtained by 
plugging (jl.9|) and (j2.11j) into the expressions ()3.8|) and (j3.11|) for the coefficients a n (y,u; A). 

THEOREM 3.4. Suppose that an electric potential V(x) and a magnetic field F(x) are 
C°° -functions and that they admit the asymptotic expansions and (jl.l^) where Vj(x) 

and Fj(x) are homogeneous functions of orders —pj and —rj = —pj — 1, respectively, where 
1 < pi < p 2 < — Let the magnetic potential A(x) be defined by equalities (\2. 6|) - [\2. &|) . 
Then, the symbol a(y,u; A) of the PDO S(X) admits the expansion into the asymptotic sum 

oo oo 

a(y,u;X) ~ 1 + Yl a ^™m,j2,-,jn(y,u; \), (3.20) 

n=l m=0 jija,— ,jn 

where jk = 1,2,... for all k = 1, . . . , n, m = 0, 1, . . the functions o J n,m;j 1 ,j 2 ,--;j n (y^ w\ A) only 
depend on , V J2 , . . . , Vj n and are homogeneous of order 

n-m- p h - p h - p jn 

in the variable y. In particular, 

a li0; j(y, u; A) = R(y, uj; A; V,-) (3.21) 
is a homogeneous function of order —pj + 1. 

COROLLARY 3.5. //, up to functions from the Schwarz class, V(x) and F(x) are homo- 
geneous (for large \x\) functions of orders —p and — r = —p — 1, then 

oo oo 

a(y,cj;\) ~ 1 + ^ 

(j/,w;A), 

n=l m=0 

where a n ^ m is of order n in V (linear, quadratic, etc.) and is a homogeneous function of 
order n — m — np of the variable y. 

4 The Inverse Problem 

4.1 Inversion of the Radon transform. Our approach to the inverse problem relies 
on the two-dimensional Radon transform (see, e.g., J!]). Let us recall it briefly here. For 
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v G S P (M. 2 ), p > 1, the Radon (or X-ray, which is the same in the dimension two) transform 
is defined by the formula 

/oo 
v(ujt + y)dt 1 u) G §, (uj,y) = 0. (4.1) 
-oo 

Obviously, r(u,y) = r(—u,y). The Fourier transform v of v and hence the function v itself 
can be recovered in the following way. Let uj^ be one of the two unit vectors such that 

/oo 
e-^ s r(^,si;v)ds, t = £\t\- 1 - ( 4 - 2 ) 
-oo 

We apply this method for the reconstruction of a homogeneous function V G C°°(lR d \{0}) 
of order — p < — 1 from its X-ray transform (jl.7|) known for all uo G S d_1 and y G IT.,, y ^ 0. 
For an arbitrary xGlS^ {0}, consider some two-dimensional plane A x orthogonal to x and, 
for y G Aa;, set v x (y) = V(x + y). Then for all u G A x , \u\ — 1, and all y G A x , (u, y) = 0, 

r(u,y;v x ) = R e (uj, x + y;V). (4.3) 

Since x + y ^ 0, the function ^ G iS _p (R 2 ) so that we can recover the function v x and, in 
particular, v x (0) = V(x) by formula (|4.2|) . This procedure is used for the reconstruction of 
the asymptotics of the electric potentials. 

In the magnetic case we are given only the integral (jl.8|) for all u G and y G n^, 
y 7^ 0. Since this integral is zero if A(x) = grad0(x) and <f)(x) — > as |x| — > oo, we 
cannot hope to recover A from this equation. Nevertheless the corresponding magnetic field 
F(x) = curlv4(x) can be recovered. Assume again that A G C°°(lR d \ {0}) is a homogeneous 
vector- valued function of order — p < —1. Let us consider one of the components of F(x), for 
example, F^ 12 \x). We will first show how F^ 12 \x) can be reconstructed everywhere except 
the plane L 12 where x 3 = . . . = Xd = 0. Let u = (001,002, 0, . . . , 0) be any unit vector in the 
plane L 12 , let v = (—oo 2 , 0J\, 0, . . . , 0) be the unit vector obtained by rotating u in the plane 
L 12 by the angle 7r/2 in the counter-clockwise sense and let x = (0, 0, £3, . . . , x n ) 7^ be an 
arbitrary vector that is orthogonal to L 12 . Set f~ (y) = F^ 12 \x + y) for y G L 12 . It is easy 
to see that 

r(uj,y; f^ 12) ) = -dR m (uj,sis + x;A)/ds, y = sv. (4.4) 
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Indeed, since is invariant under rotations in the plane L12, it suffices to check (14. 4 j) for 
the case uj% = 1, w% — when (|4.4|) reads as 

/OO £ /-oo 

F (12) (t, s, x 3 , . . . , x n )dt = -tt ^4 (1) (t, S,X 3 ,..., X n )dt. 
00 OS J — 00 

For the proof of this relation, we have only to use the definition F^ = dA^ Jdt — dA^ /ds 
and the fact that the integral of dA^ /dt is zero. Finding expression (|4.4J) . we can recover 
F^ 12 \y + x) for all y G L12 by formula (|4.2jh By virtue of the homogeneity of the function 
i 7 ^ 12 - 1 , this yields i 7 ^ 1 ^^) everywhere except the plane L12. Thus, the magnetic field F(a;) 
is reconstructed for all x such that all coordinates Xi,X2 ■ ■ ■ ,Xd are not equal to zero. This 
condition can be achieved for any arbitrary x 7^ by a suitable choice of the coordinate 
system. 

Formula (|4.4J) can of course be rewritten in the invariant way. For example, in the case 
d = 3 it means that, for arbitrary w,ng§ 2 , (u, n) = 0, 

/oo /*oo 
(n, curl A{tuo + x))<it = (c<j A n, \7 X / (a;, A(tu + x))dt), 
-00 •/ —00 

where the symbol " A " means the vector product. 

Finally, if only the combination (|1.2J) is known, then using that R e and R m are even and 
odd functions of u G § d_1 , respectively, we obtain that 

R e (y,co;V) = ik(R(y,u;\;V) + R(y,-u;\;V)), 
R m (y,u;A) = i2-\R(y,-uj;\;V)-R(y,u;\;V)). 

This allows us to reconstruct the functions V and F by the formulae given above. 
In particular, we obtain 

PROPOSITION 4.1. Let V and F be homogeneous functions on M. d ,d > 3, of orders 
—p < — 1 and —r = —p — 1 . If R(y, ui; A; V) = for all uj G S^ 1 and y G 11^, y 7^ 0, then 
V(x) = F(x) = 0. 

4.2 Reconstruction of the potential at infinity. Now we are in a position to recon- 
struct, under the assumptions of Theorem 3.4, all asymptotic coefficients Vj and Fj in (|1.9|) 
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and (jl.lUj) from the kernel s(u, u; A) of the scattering matrix known, up to a C 00 (S d_1 x § d_1 )- 
function, for some A > and all v, ui G S d_1 . Thus, we assume that V and F admit the 
asymptotic expansions (jl.9|) and (jl.lUj) . Then f)2.12j) holds but the functions Vj, and, in par- 
ticular, their orders of homogeneity —pj, are not known. Let us define the function a(y, u; A) 
by formula (|3.17j) . that is a(y,u;X) is the right symbol of the pseudodifferential operator 
S(X). This function admits an asymptotic expansion as in ()3.13|) where do = 1 and a n are 
homogeneous functions of order — p/. and < p\ < p-i < — Our goal is to solve equation 
(I3.18|) . where the function a is known with respect to V. Of course both sides of (j3.18j) are 
defined up to functions from the Schwarz class S. Under the assumptions of Theorem 3.4 
the orders of homogeneity in the left- and right-hand sides of this equation, as well as the 
terms of the same order, should coincide. 

Let us compare terms of the highest order in (j3.18|) . Taking into account relation (j3.16j) . 
we see that p\ = p\ + 1 and 

R(y,u;\;V 1 ) = a 1 (y,w',\). (4.5) 

This allows us to reconstruct the coefficients V\ of the electric potential and F\ of the 
magnetic field by the formulae of the previous subsection. Then we find A\ by formulae 

- (Q. 

Below it is convenient to introduce the following notation. Suppose that some function 
/ admits the expansion in the asymptotic series (|2.9|) of homogeneous functions. By taking 
sums of such functions it cannot be excluded that some terms will be equal to zero. Therefore 
we define /" as the highest order homogeneous term /& in ()2.9|) that is not identically zero. 
For example, /" = if / G <S. 

Suppose now that we have found the coefficients for all k — 1, . . . , n — 1, n > 2. Let us 
reconstruct V n . We apply Proposition 3.3 to the functions V^ 1 -* = Yl^Zi Vj, V 1 - 2 - 1 = V — V'- 1 -' 
where p^ = p\, = p n . This yields 

n-l 

Q(v)-g(^v J )G5-^ +1 - e , e>o, 

3=1 

and thus, for an arbitrary p n , this term can be neglected compared to R(V n ). All terms 
R(Vj), j > n + 1, are also negligible compared to R(V n ). Therefore rewriting equation 
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(irrsi) as 



n-1 



n-1 



i?(v n )+ R(y j ) + (Q(v)-Q(J2 v ^ 



a-l-T(£,Vj) 




and selecting terms of the highest order, we obtain that 



n-1 



R(V n ) = (a-l-TC£v j )) t - 



(4.6) 



Having found R(V n ), we can, similarly to f)4.5|) . recover the functions V n , F n and, conse- 
quently, using formulae ()2.6|) - ()2.8|) the corresponding A n . 

Thus, we have arrived at our main result. 

THEOREM 4.2. Suppose that an electric potential V {x) and a magnetic field F(x) are ex- 
junctions on M. d , d > 3, and that they admit the asymptotic expansions (j i.ffi and (\l.l(Jty where 
Vj(x) and Fj(x) are homogeneous functions of orders —pj and —rj = —pj — 1, respectively, 
where 1 < p\ < p% < ■ • -. Let the magnetic potential A{x) be defined by equalities i\2. 6|) - &|) . 
Then the scattering data consisting of the kernel s(u, to'] A) of the scattering matrix at a fixed 
positive energy in a neighbourhood of the diagonal v = uj uniquely determines each one of 
the Vj(x) and the Fj(x). Moreover, the functions V\(x) and Fi(x) can be reconstructed from 
formula ( |^.<5| ) and the functions Vj(x) and Fj(x) for j > 2 can be recursively reconstructed 
from formula Q^.6| ) . 

COROLLARY 4.3. // kernel of the operator S(X) — I belongs, for some X > 0, to the 
space C 00 ^- 1 x S^ 1 ), then V G S(R d ) and A G S(R d ). More generally, let V (1 \F^ 
and V^ 2 \ F^ satisfy the assumptions of Theorem 4.2. Suppose that s^ l \X) — s^(A) G 
C7 00 (S rf - 1 x S '- 1 ) for some A > 0. Then V (1) - V (2) and F« - F (2) belong to the Schwarz 
class S(R d ). 

Suppose that, for some n > 2, 



n-1 



a — 



i-t^v^gsou 



for all uj G S^" 1 . Then it follows from fO| that Vj = Fj = for all j >n. 
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We emphasize that, unlike p\ = p,\ + 1, other orders — p n are not determined by the 
sequence pi, f-iz, ■ ■ ■ only, in particular, because a cancellation of different terms in the right- 
hand side of ()4.6|) might occur. For example, in 

R(y 2 ) = (a-l- T(Vi))« = (a - 1 - ax - Q(Vi))« (4.7) 

a cancellation of terms in a — 1 — ax and Q(Vi) cannot be excluded. Nevertheless — p 2 should 
be smaller than the maximal order of these two series. The order of a — 1 — a\ equals — p 2 
and the order of Q(Vi) equals or is smaller than — min{pi, 2pi — 2}. This yields the bound 

p 2 > min{/! 2 , pi, 2pi - 2} = min{p 2 , Hi + 1, 2pi}. 



4.3 A counter-example. We give now a counter-example that shows that in the 
dimension two the argument above does not work. Let us first calculate the Radon transform 
of the function V(x) = v(x)\x\~ 2 where x = G S. Making the change of variables 

t = \y\r, t = tans, we see that 

r u/ , , . f°° fy + Tu\ dr . r /2 , A . 

/ + = y / v \ —=== \ —— = \y\ I v(cossy + sm sujds 

J -co J-oo \Vt 2 + 1Jt 2 + 1 J_ n/2 

= \y\- 1 [ v(9)d9, (u,y) = 0, y + 0, if y G S ± (^), (4.8) 

JS±(<j) 

where S + (cj) (§ _ (o;)) is the half-circle passed from the point —uo to to in the counter-clockwise 
(clockwise) direction. 

It follows that the Radon transform of V is equal to zero for all lj G S, (uj, y) — 0, y ^ 0, 
if and only if 

v(#)cifl = 



'S+(w) 

for all G S. It is easy to see that this condition is equivalent to two conditions 



v (uj) = v (- u ) and / v(0)d0 = 0. (4.9) 
For example, they are satisfied if 

v(lu) = 2{u>, ujq) 2 — 1 

where ujq G S is some fixed vector. 
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COUNTER-EXAMPLE 4.4. Let d = 2, let A = and let the electric potential V be a 
C°°-function such that V(x) = v (x)\x\~ 2 for sufficiently large \x\. Assume that the function 
v satisfies conditions (|4.9|) but is not identically zero. Then the Radon transform ()4.8|) of 
v(x)|x|~ 2 equals zero for all uj G S, (u,y) = 0, y ^ 0, although v does not. According to 
Theorem 3.2, in this case 5(A) — I is the PDO from the class S~ 2 . 

We emphasize that in this example the leading singularity of the scattering amplitude 
disappears for all energies A > 0. The examples of pi] are much stronger in the sense that 
the scattering amplitude (not only its leading singularity) equals zero. On the other hand, 
this is true for one energy only. It is interesting that both the examples above and of ^T] 
are two-dimensional and the potentials decay at infinity as homogeneous functions of order 
-2. 

4.4 Uniqueness result. In this subsection we assume that the magnetic field is zero and 
that d > 3. According to Corollary 4.3 the electric potential is determined by the scattering 
matrix 5(A) for some A > up to a term from the Schwarz class. Let us show now that 
under stronger a priori assumptions this remainder can be removed. To that end we need 
the following result which is a particular case of Theorem 1 of [34J. 

THEOREM 4.5. Let the potentials V {1) and V {2) on R d , d> 3, satisfy estimates (J23) with 
p > d and let V^ 1 '^) = V^ 2 '(x) for sufficiently large \x\. Suppose that the corresponding 
scattering matrices S^\\) = (A) for some A > 0. Then V^(x) = V^ 2 \x) for all x G M d . 

Combining Theorems 4.2 and 4.5, we obtain a new uniqueness result. 

THEOREM 4.6. Let C°° -potentials V {k) , on R d , d > 3, k = 1,2, equal (finite) sums of 
homogeneous functions of order —pf^ for sufficiently large \x\. Let pj > d for k = 1, 2 and 
all j. Suppose that the corresponding scattering matrices S^(X) = S^(\) for some A > 0. 
Then VW(x) = V®(x) for all x G R d . 

Indeed, it follows from Corollary 4.3 that V^'(x) = V^ 2 '(x) up to a function from the 
Schwarz class. By our a priori assumption this function has compact support. Then it 
follows from Theorem 4.5 that, actually, this function is zero. 
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